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ABSTRACT 

We present adaptive optics (AO) near-infrared observations of the core of the Tr 14 cluster in the Carina region obtained with the 
ESO multi-conjugate AO demonstrator, MAD. Our campaign yields AO-corrected observations with an image quality of about 0.2" 
across the 2' field of view, which is the widest AO mosaic ever obtained. We detected almost 2000 sources spanning a dynamic range 
of 10 mag. The pre-main sequence (PMS) locus in the colour-magnitude diagram is well reproduced by Palla & Stabler isochrones 
with an age of 3 to 5 x 10^^ yr, confirming the very young age of the cluster. We derive a very high (deprojected) central density 
«o ~ 4.5(±0.5) X 10'* pc"^ and estimate the total mass of the cluster to be about ~ 4.3^j j x 10' Mq, although contamination of 
the field of view might have a significant impact on the derived mass. We show that the pairing process is largely dominated by 
chance alignment so that physical pairs are difficult to disentangle from spurious ones based on our single epoch observation. Yet, we 
identify 150 likely bound pairs, 30% of these with a separation smaller than 0.5" (~1300AU). We further show that at the 2cr level 
massive stars have more companions than lower-mass stars and that those companions are respectively brighter on average, thus more 
massive. Finally, we find some hints of mass segregation for stars heavier than about 10 Mq. If confirmed, the observed degree of 
mass segregation could be explained by dynamical evolution, despite the young age of the cluster. 

Key words. Instrumentation; adaptive optics - Stars: early-type - Stars: pre-main sequence - binaries: visual - open clusters and 
associations: individual: Tr 14 



1. Introduction 

Massive stars do not form in isolation. They are bom and, for 
most of them, are living in OB associations and young clusters 
( |Mafz-Apellaniz et al. 2004j). Indeed, most of the field cases 
are runaway objects that can be traced back to their natal clus- 
ter/association ( |de Wit et al.||2005] l. Even the best cases of field 
massive stars are now questioned in favour of an ejection sce- 
nario ( [Gvaramadze & Bomans 2008 ). 

One of the most striking and important properties of high- 
mass stars is their high degree of multiplicity. Yet accurate ob- 
servational constraints of the multiplicity properties and of the 
underlying parameter distributions are still lacking. These quan- 
tities are however critical as they trace the final products of high- 
mass star formation and early dynamical evolution. In nearby 
open clusters, the minimal spectroscopic binary (SB) fraction 
is in the range of 40% to 60% ( Sana et al.| [Sodg] [2009] [2010) 1, 
which is similar to the 57% SB fraction observed for the galactic 
O-star population as a whole (Mason et al. 2009 1. 

While spectroscopy is suitable to detect the short- and 
intermediate-period binaries (P < 10 yr), adaptive optics (AO) 
observations can tackle the problem from the other side of the 



separ ation range (see e.g. discussion in [Sana & Le Bouquin 



|2009 1. As an example. Turner et al. ( 2008| l obtained a minima 



* Based on observations obtained with the MCAO Demonstrator 
(MAD) at the VLT Melipal Nasmyth focus (ESO public data re- 
lease). Tables |3] and [5] are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
[http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-H A/| 



fraction of massive stars with companions of 37% within an an- 
gular separation of 0.2 to 6". However, their survey is limited to 
objects with declination 6 > -42°. It is thus missing some of the 
most interesting star formation regions of the Galaxy, like the 
Carina nebula region. 

In this context, we undertook a multi-band NIR AO 
campaign on the main Carina region clusters with the Multi- 
Conjugate Adaptive Optics (MCAO) Demonstrator (MAD, 
Marchett Fet al.[ \2001} . Beside deep NIR photometry of the 
individual clusters, our survey was designed to provide us with 
high-resolution imaging of the close environment of a sample 
of 60 O/WR massive stars in the Carina region. Unfortunately, 
the bad weather at the end of the second MAD demonstration 
run in January 2008 prevented the completion of the project. 
Valuable H and K5 photometry of the sole Tr 14 cluster could 
be obtained. The 2' field of view (fov) still provides us with 
high-quality information of the surrounding of ~30 early-type 
stars with masses above 10 Mq. It also constitutes the most 
extended AO mosaic ever acquired. 

Located inside Carina at a distance of 1.5-3.0 kpc, Tr 14 is 
an ideal target to search for multiplicity around massive stars be- 
cause it contains more than 10 O-type stars and several hundreds 
of B-type stars ( Vazquez et al.[ [l996). Large differences in the 
distance to Trumpler 14 arise from adopting different extinction 
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Table 1. Field centering (EC.) for on-object (Tr 14) and on-sky 
observations and coordinates of the natural guide stars (NGSs). 







RA 


DEC 


V mag 






Tr 14 F.C. 


10:43:55.00 


-59:33:03.0 








Sky F.C. 


10:43:06.53 


-59:37:46.1 








NGSl 


10:43:59.92 


-59:32:25.4 


9.3 






NGS2 


10:43:57.69 


-59:33:39.2 


11.2 






NGS3 


10:43:48.82 


-59:33:24.8 


10.7 




Table 2. Log of the MAD observations of Tr 14. 


DP 


RA 


DEC 


DIT NDIT 


Mma 


Tot. 






Trumpler 14 observations 






#1 


10:43:56.0 


-59:32:46 


2s 30 


2x14 


28 min 


#2 


10:43:56.5 


-59:33:29 


2s 15 


2x8 


8 min 


#3 


10:43:53.5 


-59:33:28 


2s 15 


2x8 


8 min 


#4 


10:43:53.5 


-59:32:41 


2s 15 


2x8 


8 min 




Sky field observations (MCAO in open loop) 




#1 


10:43:07.5 


-59:37:29 


2s 30 


8 


8 min 


#2 


10:43:08.0 


-59:38:12 


2s 15 


8 


4 min 


#3 


10:43:05.0 


-59:38:11 


2s 15 


8 


4 min 


#4 


10:43:05.0 


-59:37:24 


2s 15 


8 


4 min 



DP #4 



NGSl 




Fig.l. 2MASS K band image of Tr 14. The cross and the 
large circle indicate the centre and the size of the 2' MAD 
field of view. The four 57"x57" boxes show the position of the 
CAMCAO camera in the adopted 4-point dither pattern, while 
the selected NGSs are identified by the smaller circles. 




Fig. 2. Averaged FWHM and Strehl ratio maps as computed over 
the MAD field of view (yellow circle). The three red crosses 
show the locations of the NGSs. North is to the top and East to 
the right. 



laws and evolutionary tracks (Carraro et al. 2004[ l. Differences 
in distance can partially account for diverse estimates of mass 
and structural parameters. Its mass was first estimated to be 
2000 Mq (|Vazquez et al. 19961. However, the photometry used 
by these authors barely reached the turn-on point of the pre-main 
sequence (PMS), while they extrapolated the mass assuming a 
Salpeter initial mass function (IMF). More recently, Ascenso 
|et al.| ( |2007| l used much deeper IR photometry, which revealed 
the very rich PMS population and provided a more robust mass 
estimate of 9000 M q. [Vazquez et alT ( 1996 1 reported a core ra 



and detected for the first time a core-halo structure, which is typ- 
ical of these young clusters (e.g., Baume et al. 2004[ l. Tr 14 is 



dius of 4.2 pc, while [Ascenso et al.|p007| l revised it to 1.14 pc. 



indeed very young, not yet relaxed and has been forming stars in 
the last 4 Myr ([V azquez etaLl[1996 |l. 

The layout of the paper is as follows. Sections |2] and [3] de- 
scribe the observations, data reduction and photometric analy- 
sis. Section |4]presents the NIR properties of Tr 14 and discusses 
the cluster structure. Section|5]analyses the pairing properties in 
Tr 14. Section |6]describes an artificial star experiment designed 
to quantify the detection biases in the vicinity of the bright stars. 
It also presents two simple models that generalise the results of 
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Fig. 3. False colour image of the 2' fov of the MAD observations (blue is H band; red, /Ts band). North is to the top and East to the 
left. 



the artificial star experiment. As such, it provides support to the 
results of this paper. Finally, Sect.|7]investigates the cluster mass 
segregation status, and Sect.|8]summarizes our results. 

2. Observations and data reduction 

The MAD instrument is an adaptive optics facility aiming at 
correcting for the atmospheric turbulence over a wide field-of- 
view, and as such constitutes a pathfinder experiment for MCAO 
techniques. Briefly, MAD relies on three natural guide stars 
(NGSs) to improve the image quality (IQ) over a 2'-diameter 
fov. Optimal correction is reached within the triangle formed by 
the three NGSs although some decent correction is still attained 
outside, mostly depending on the observing conditions and on 
the coherence time of the atmospheric turbulence. 

The CAMCAO IR camera images a 57"x57" region in the 
MAD fov and is mounted on a scanning table, so that the full 2'- 
diameter fov can be covered with a 4- or 5 -point dither pattern. 
The detector used is an Hawaii 2kx2k, yielding an effective pixel 
size on sky of 0.028". 

Because of the constraints imposed by the geometry and the 
magnitudes of the NGSs as well as by the brightness limit of the 
detector (typically > 8), the Carina clusters turned out to be 



ideal targets for our purposes. Combined with the large collect- 
ing area of an 8-m class telescope, MAD was offering a unique 
opportunity to collect the missing high-spatial resolution obser- 
vations to characterize a statistically significant set of massive 
stars. 

On the night of January 10, 2008 during the second Science 
Demonstration (SD) run, the MAD team acquired H and Ks 
band observations of Tr 14. Because of the mentioned con- 
straints on the choice of the NGSs, the MAD pointing was offset 
from the cluster centre by about 0.4' to the W-SW and a 4-point 
dither pattern was used to cover the (almost) full 2' diameter 
fov (Fig. [TJ. We obtained 28 images for a total exposure time 
of 28 min on the central field of the cluster (DP#1) and eight 
30sec images in the three remaining DPs. The size of the jitter 
box was 10". We further followed a standard object-sky-object 
strategy. The jitter and dither pattern for the on- and off-target 
observations were identical, although the sky observations were 
obtained without AO correction. The sky field, located about 5' 
SW from the cluster centre, is one of the few IR-source depleted 
regions of the neighbourhood. The journal of the on- and off- 
target observations is summarized in Table |2] For each dither 
position (Col. 1), Cols. 2 and 3 list the coordinates of the four- 
point dither pattern. The detector integration time (DIT) and the 
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Fig. 4. Close-up view of the Ks band image around selected mas- 
sive stars. The image size is 200x200 pixels, corresponding to 
approximately 5.5x5.5" on the sky. 



number of repetitions at each jitter position (NDIT) are given in 
Cols. 4 and 5. Finally, Cols. 6 and 7 respectively indicate the 
number of images and the total integration time spent on each 
dither position. 

The data were reduced with the iraf package. All science and 
sky frames were dark-subtracted and flatfielded with the calibra- 
tions obtained by the SD-team in January 2008. A master-sky 
was created for each dither position (DP) by taking the median 
of the corresponding sky images. The few visible sources in the 
sky images were manually masked out before computing the me- 
dian sky. This master-sky was subtracted from the individual sci- 
ence images. While the point spread function (PSF) photometry 
(see Sect. [3]l was performed on the individual images, we also 
combined the images into a 2' diameter mosaic that was used 
to estimate the overall IQ of the campaign. We note the pres- 
ence of electronic ghosts offset by « x 128 pix from a bright 
source (with n - ±1, ±2, +3, ...) along the detector reading di- 
rection. We thus systematically masked out the corresponding 
sub-regions. Thanks to the jittering and to the fact that two con- 
tiguous CAMCAO quadrants do not have the same read-out ori- 
entation, most of the masked zones were recovered when com- 
bining the different frames. 

Ambient conditions during our observations were as fol- 
lows. The R band seeing was varying between 0.9" and 1.8", 
corresponding to a band seeing between 0.7" and 1.5". The 
coherence time was in the range of 2 to 3 ms. Even though the 
ambient conditions were clearly below average for the Paranal 
site, the MCAO still provided a decent improvement with an 
full-width half maximum (FWHM) of the PSF of 0.2" over 
most of the 2' fov. The corresponding Strehl ratio was estimated 
in the range of 5-10% (Fig.|2]i. Figure [3] displays a false colour 
montage of our data on Tr 14, while Fig. |4]presents a close-up 
view on a 5x5" region, with the aim to emphasize the shape and 
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Fig. 5. Comparison of the H (upper panel) and Ks band (lower 
panel) photometry obtained in this paper (MAD) with that of 
[Ascenso et al. (2007 SOFI). The squares indicate the stars used 
to compute the zero point difference. 
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Fig. 6. H and Ks band l-cr error bars on the PSF photometry 
plotted versus the star magnitude. 



smoothness of the PSF, even on the co-added images. 

3. PSF photometry 

Stellar photometry was obtained with the PSF fitting technique 



using the well- tested daophot/allstar/allframe (Stetson 1987 
1 199 4) packages. The advantage of using allframe is that it em- 
ploys PSF photometry on the individual images, thereby it ac- 
counts better for the varying near infrared sky and seeing con- 
ditions. The latter have indeed a significant impact on the qual- 
ity of the AO correction and thus of the actual IQ of the data. 
The PSF of each individual image was generated with a PENNY 
function that had a quadratic dependence on position in the 
frame, using a selected list of well isolated stars. 
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ID RA 


DEC. 


H 


cth 




0"Ks 


HH:MM:SS.sss 


±DD:AM:AS.s.s 











1 


10:43 


46 


953 


-59 


33 


18 


40 


12 


192 


0.009 


11.420 


0.008 


2 


10:43 


47 


013 


-59 


32 


43 


47 


16 


406 


0.011 


15.794 


0.017 


3 


10:43 


47 


017 


-59 


32 


51 


06 


15 


233 


0.008 


14.668 


0.010 


4 


10:43 


47 


017 


-59 


33 


09 


70 


11 


828 


0.010 


10.738 


0.009 


5 


10:43 


47 


026 


-59 


32 


42 


31 


11 


506 


0.008 


11.301 


0.008 


6 


10:43 


47 


090 


-59 


33 


24 


45 


17 


031 


0.032 


16.428 


0.011 


7 


10:43 


47 


095 


-59 


33 


17 


39 


15 


317 


0.012 


14.687 


0.008 


8 


10:43 


47 


125 


-59 


33 


12 


57 


17 


673 


0.036 


16.884 


0.015 


9 


10:43 


47 


186 


-59 


33 


21 


88 


14 


770 


0.010 


13.982 


0.007 


10 


10:43 


47 


241 


-59 


33 


17 


26 


15 


511 


0.014 


15.025 


0.009 


11 


10:43 


47 


309 


-59 


33 


06 


09 


13 


557 


0.012 


12.963 


0.010 


12 


10:43 


47 


382 


-59 


33 


41 


01 


14 


883 


0.007 


14.508 


0.006 


13 


10:43 


47 


430 


-59 


33 


39 


45 


16 


863 


0.019 


16.231 


0.008 


14 


10:43 


47 


464 


-59 


33 


21 


92 


14 


537 


0.012 


14.070 


0.008 


15 


10:43 


47 


497 


-59 


32 


40 


79 


14 


600 


0.008 


13.979 


0.009 


16 


10:43 


47 


505 


-59 


33 


39 


95 


14 


184 


0.006 


13.490 


0.006 


17 


10:43 


47 


516 


-59 


33 


10 


58 


14 


125 


0.012 


13.773 


0.009 


18 


10:43 


47 


519 


-59 


33 


26 


10 


17 


495 


0.036 


16.973 


0.013 


19 


10:43 


47 


534 


-59 


33 


03 


60 


15 


146 


0.008 


14.545 


0.006 


20 


10:43 


47 


539 


-59 


33 


18 


65 


18 


158 


0.043 


17.211 


0.019 



note: The full version of the table is available in the electronic version of the paper or through the CDS: |http://cds.u-strasbg.&7| . 



The calibration of the instrumental H and MAD data was 
done by direct comparison with the JHKs calibrated SOFl cat- 
alogue of |Ascenso et al.| ( |2007| l. Of particular importance is the 
absence of (i) any colour term between the two photometric sys- 
tems over a six-magnitude range (Fig.|5]l, and (ii) spatial system- 
atics between the 4-single MAD pointings and the mean zero- 
point difference with respect to the SOFl data. Open squares in 
Fig. |5] highlight the stars used to estimate the mean offset be- 
tween the two systems, which were selected by applying a So- 
clipping around the mean zero-point difference. We obtained 



A/Ts = -0.008 + 0.122 
AH ^ -0.007 ±0.101 



(1) 
(2) 



Taking into account that some objects might be variable and 
that some others were likely not resolved by SOFl, we conclude 
that there is an almost perfect agreement between the two sets 
of measurements. We therefore continue without applying any 
correction to our photometry. After manually cleaning the hand- 
ful of double entries, our photometric catalogue contains 1955 
stars, most of them brighter than - 19 or // = 20 mag. A 
sample of the catalogue is given in Table |3] where Col. 1 in- 
dicates our internal identifier. Cols. 2 and 3 give the equatorial 
coordinates and Cols. 4 to 7 provide the H and Ks band magni- 
tudes and the Icr error bars. An overview of the accuracy of the 
magnitude measurements according to the object brightness is 
provided in Fig.|6]. The uncertainties are typically of cr a; 0.005 
for stars brighter than = 15 mag (H = 17 mag resp.) and 
reach cr ^ 0.02 at Ks - 18 mag {H - 19 mag resp.). 



diagram (CMD) of Tr 14 for dither positions DP #1 and 4, and 
for the full fov (DP#l-4) and compares it with the main sequence 
(MS) and PMS locations given the adopted cluster distance and 
reddening. Figure|7]further provides an overview of the light-to- 
mass conversion scale used in the rest of this paper 

While most of the stars brighter than Ks - 12 agree well with 
the MS of Lejeune & Schaerer ( 2001| l, the vast majority of the 
fainter stars (Ks > 14) are still m the PMS stage. A comparison 



with the PMS isochrones of Palla & Stabler (1993i suggests a 



contraction age younger than 1 Myr, and possibly as young as 
0.3-0.5 Myr. At this age, the transition between the PMS and the 
MS occurs for stars with masses in the range of 4 Mq to 8 Mq. 
Although the PMS isochrones are still affected by uncertainties 
in the colour transformation, our data clearly suggest that the 
core of Tr 14 has undergone a very recent starburst event, during 
which most of the low- and intermediate-mass stars in Tr 14 have 
been formed. Using a larger distance and/or a larger reddening 
(e.g., Ascenso et al. 2007| l would result in even younger ages, 
which is the reason why we conservatively decided to adopt the 
formerlCarraro et aljresults. 

Without control field observations, we cannot quantitatively 
estimate the contamination of the CMD by field stars. Yet, the 
CMD from DP#4 shows the least structure and provides some 
qualitative estimate of the maximum contamination suffered by 
the central part of the cluster. It indicates that the dispersion 
around the 3 x 10^ yr isochrone is mostly real. Comparing them 
with PMS isochrones of different ages, we estimated the dura- 
tion of the starburst to be a couple of 10^ yr at most. 



4. A MAD view of Tr 1 4 

4.1. Colour-magnitude diagram 

Because we only have H and Ks band observations, our data 
alone cannot altogether constrain the reddening, the distance and 
the age of the stellar population. Below, we adopt the results of 
Carraro et aL| ( [2004| ): d ^ 2.5 + 0.3 kpc, DM = 12.3 + 0.2 
and Ay - 2.0 ± 0.13. Figure |7] presents the colour magnitude 



4.2. Foreground/background contamination 

Highly reddened stars seen in the CMD cannot be reproduced 
by any evolutionary track and are thus deeply embedded clus- 
ter members, background objects whose reddened colour results 
from ISM absorption or foreground cool stars. Similarly, very 
blue stars in the CMD are unlikely members of Tr 14. 

Before statistically estimating the likelihood of a given vi- 
sual pair to be physically bound (Sect. |5]), we first used the cur- 
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Fig. 7. L eft: Tr 14 CMP for DP #1. M iddle: CMD f or DP #4. Right: c omplete Tr 14 CMD (DP#1 to 4). The dashed lines show the 
MS from Lejeune & Schaerer (2001 1 and from Palla & Stahler ( 1993 1 for stars with masses M > 6 Mq and M < 6 Mq respectively. 



Diamonds and left-hand labels indicate the MS masses. The plain lines show the log{age/yr) -5.5 PMS isochrone from Palla & 



Stahler for PMS stars with masses between 0. 1 and 6.0 Mq. Squares and right-hand labels indicate the corresponding PMS masses. 
12.3 and Ay = 2.0 were adopted as discussed in Sect. 



DM 



4.1 



10 



20 





























■ + 






















H-Ks 



Fig. 8. Tr 14 CMD. The plain lines delimit the adopted locus of 
cluster members (see text) while the stars show the objects not 
considered in Sect. |5] 



rent results to clean our catalogue from improbable cluster mem- 
bers. For bright stars {K^ > 12 mag), we required the H - 
colour to be within 0.2 mag from the expected locus of MS stars 
with M > 6 Mq. For stars fainter than Ks = 12 mag, we adopted 
the following relations as the left and right H - Ks limits in the 
CMD (Fig.[|): 



(H-Ks)Mt = 0.043(/:s - 12), 
(H - ^:s)nght = 0.3 - 2.0/(Ks - 20). 



(3) 
(4) 



Both relations are obtained as an approximate to the +2cr inter- 
val around the typical H - Ks colour of an 0.3 Myr PMS star 
for a given Ks magnitude. Applying those criteria, we excluded 
12 bright and 233 faint stars. This represents about 12.5% of 
our total sample. Below, we also focus on stars brighter than 
Ks = 18 mag. This excludes an additional 198 stars (a!lO% 



o + 
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Fig. 9. Distribution of the reddest and bluest stars in Trl4. The 
dot and the open circles show the location of the faint {Yl <Ks < 
18 mag) and bright {Ks > 12 mag) stars respectively. Asterisks 
and crosses identify the bluest and reddest faint stars, while the 
filled circles trace the reddest bright stars. For the latter, the sym- 
bol size is proportional to the Ks magnitude. The dashed circle 
shows a r diameter region centered on HD93129A. North is at 
the top and East to the left. 



of our sample), resulting in a list of 1495 likely members with 
/Ts < 18 mag. 

Figure |9] displays the spatial distribution of the rejected blue 
and red stars. The faint reddest stars agree well with a random 
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4.3. Cluster structure 
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Fig. 10. Tr 14 surface number density distributions for the full 
cluster population (circles) and for the PMS (triangles) and MS 
(diamonds) populations. Upper and lower panels show the den- 
sity profiles obtained, respectively, before and after applying the 
colour criteria of Fig. [8] In both cases only stars with magnitude 
/Ts < 18 mag have been considered. Different lines show the 
best-fit EFF87 profiles described in Table |4] 



spread in the field. The faint bluest ones are mostly located in 
the North and East edges and suggest larger color uncertainties 
in that zone. The surface density of the brighter red stars how- 
ever shows a clear enhancement that correlates with the central 
part of the cluster. This suggests that the corresponding stars 
are rather highly reddened objects associated with Tr 14. We 
acknowledge that we might thus have rejected some members 
of Tr 14. Preserving the homogeneity of the reddening proper- 
ties of the bright star sample is however more important and we 
argue that the resulting few extra rejections will not aff'ect the 
statistical companionship properties discussed in Sect. |5] 



Adopting the cluster centre as defined by Ascenso et al. ( 2007| l, 
we computed the radial profile of the star surface density 
(Fig. 10 1. Ascenso et al. (2007i proposed a core-halo structure 



with a respective approximate radius of 1' and 5'. Although our 
data set covers a limited area, there is no indication of a transition 
between the two regimes. We can thus consider that our fov is 
strictly dominated by the core of the cluster. The cluster param- 
eters derived below thus only apply to the core of the core-halo 
structure. 

To better quantify the surface density variations, we fitted 
|Elson et al.| ([r987 hereafter EFF87) profiles, better suited for 
young open clusters than King profiles, to all three populations. 
Following EFF87, we adopt the notation 



S(r) = Sod + r'/a') 



(5) 



where Sq is the surface number density at the centre. With the 
best-fit parameters (Table |4| we also computed the radius, r^, 
where the surface number density drops to half its centre valuq^ 
and the (deprojected) central number density, n^, following Eqs. 
22 and 13b of EFF87, respectively. The central density obtained, 
«o - 4.5;^U 5 X 10"* pc"-', is very high. 

Integrating the surface number density profile to infinity and 
assuming an average stellar mass of 0.64 M© as suitable for a 
|Kroupa| ( |2001) ) IMF, our results indicate an asymptotic mass of 
~ 10 Mq, with 30% of the mass within the inner parsec 

of the profile. This value should be taken as an upper limit to the 
mass in the cluster core because of the possible contamination 
of the fov by field stars. Our mass estima te is somewhat lower 
than the value of 9 x 10^ Mq obtained by | Ascenso et al. (2007 1 
based on mass-function considerations. Yet the value of lAscensol 



et al. I falls within the upper limit of our uncertainties after cor- 
recting for the different assumption on the distances. Part of the 
difference in the mass determination might also be caused by our 
profile only probing the core region of the core-halo structure, 
which likely has a steeper profile than the halo. 

As a first attempt to search for differences between the low- 
and high-mass star properties, we also computed the density 
profiles of two sub-populations in Tr 14: the MS stars {K^ < 
12 mag) and the PMS stai's (12 < K5 < U mag). From Fig.[TO| 
(upper panel), the PMS stars occupy a larger core but display 
a steeper decrease than the MS population. Because the pop- 
ulation of the cluster is dominated by low mass stars, there is 
little difference between the number density profile of the PMS 
stars and that of the whole cluster. The more massive MS stars 
{K5 < 12 mag) however seem to be slightly more concentrated 
towards the cluster centre (see also discussion in Sect. [7]). As 
mentioned earlier, the PMS profile seems to display a steeper 
slope, but Table |4]reveals that this difference is not very signifi- 
cant (only at the 1 .80" level). 

As a second step, we also re-computed the density profiles 
after applying the colour and magnitude criteria defined in the 
previous section (Sect. 4.2 1, thus focussing on the most proba- 



ble members. The cluster core radius is found to be larger and 
the profiles display a significantly steeper decrease with radius 
(Fig.[TOj lower panel). The central surface density is also reduced 
by a factor 2.6 and the deprojected central number density by a 
factor 4.7 (Table|4]i. Following these adjustments, the asymptotic 
mass of the core is reduced to 1.4^g3 x lO^* Mq. Because of the 



' This value is often referred to as the core radius. Because Tr 14 dis- 
played a core-halo structure, should be understood as the core radius 
of the core itself. 
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Table 4. Best-fit EFF87 parameters for different populations in Tr 14. 



Population A^s range Eg [pc "] a [pc] ~y [pc] ;io [pc 









Without colour selection 










MS 
PMS 
All 


6- 
12- 
6- 


12 
-18 
18 


1077 ± 567 
10330 ±770 
11116 + 817 


0.05 + 0.03 
0.13 ±0.01 
0.13 ± 0.01 


1.66 ± 

2.17 + 

2.18 ± 


0.27 
0.08 
0.08 


"■""-0.03 

1 1 

"■^ ^-0.01 

11+""" 


0071+10497 
^"'^-6376 

40715+501'' 
44665+'^''° 








With colour selection 










MS 


6- 


12 


3097 ± 14013 


0.01 ± 0.03 


1.25 ± 


0.39 


"■"^-0.05 




PMS 
All 


12- 

6- 


-18 
18 


3977 ± 289 
4321 ± 316 


0.32 ± 0.04 
0.30 ± 0.03 


3.26 ± 
3.10± 


0.19 
0.17 


"■^-'-0.03 
"■^^-0.03 


8348![57Q 

Q4Q1+1394 



note: The uncertainties on Eq, a and 7 are the Icr error bars. The uncertainties and iiq were computed using Monte Carlo simulations assuming 
a Gaussian distribution of the errors on So, a and y. The quoted values give the 0.68 confidence interval. 
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Fig. 11. Cumulative distribution of the distance to the closest 
neighbour (limin)- The horizontal dotted line shows the size of 
our sample. 

sharper slope of the profile, one now finds 75% of the total mass 
in the inner parsec. Interestingly, the more massive stars show 
no core-structure, and their density profile is well represented 
by a simple power-law. This is somewhat comparable to what 
Campbell et al. (2010l found for the massive stars in R136. 



5. Companion analysis 

5.1. General properties 

With about 1500 likely members in a 2' -diameter fov, the mean 
surface density is 477 src/arcmin^, or 0.133 src/arcsec^. The 
closest pair detected in our PSF photometry is separated by 
0.24" and half the sources have a neighbour at no more than 
1.25" (Fig 111. Figure 12 illustrates empirically the maximum 
reachable flux contrast between two close sources as a func- 
tion of their separation. The magnitude differenc^A/Ts roughly 
scales as the cubic root of the separation: 



AKs = 6(d - 0.24) 



1/3 



(6) 



Below, the brightest star of a pair is adopted as the primary. We 



adopt the convention AKs = ^s.scc - ^■Ss.prim ^ 0, where A's.prim 
A's.scc are the primary and secondary magnitudes respectively. 



and 



Fig. 12. Magnitude difference AK^ of detected visual pairs as a 
function of the distance d. The plain line is defined by Eq.|6] 

Flux contrasts of AKs = 1/3/5/8 are reached at separations d = 
0.25/0.5/1.0/2.0" respectively. 

5.2. Chance alignment 

Because of the high source surface density, the number of pairs 
quickly increases with separation. To quantify the chance that 
an observed pair results from spurious alignment, we followed 
the approach of Duchene et al. (2001) . We define Pbound as 
the complementary probability to the one that a given pair 
(^s.prim, ^s.sec), with a Separation d, occurs by chance: 



^'bound = exp 



- z 

Ks=Ks,pL-jii 



(7) 



where R is the field radius, taken as 60". is the actual surface 
area of a star of magnitude Kg, and is given by 



WK,^n{d'-dl„{Ks-Ks,,nJi).. 



(8) 



where dmin is the minimal separation at which a star of magni- 
tude K5 can be detected in the neighbourhood of a star of mag- 
nitude /Tg prin,. It is estimated by inverting Eq.|6]thus 



fl'™„(A/:s) = (|A/:sl/6)'+0.24. 



(9) 
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Fig. 13. Probability that a pair is physically bound Pbound plotted 
vs. the binary separation d. The plain curve delimits the locus 
were no pairs are found and corresponds to Pbound = exp(-(t/ - 
0.24)2/2). 
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Fig. 14. Same as Fig. 12 where the population of likely bound 
systems has been over-plotted with red (inbound > 0.90) or blue 
squares (/"bound > 0.99). Plain red and blue lines identified the 
typical locus of the two populations. 



As expected, Fig 
physically bound pairs 



13 shows that the likelihood of finding 
/"bound, decreases steeply with separa- 
tion. Typically, only pairs with a brightness ratio close to unity 
or with a separation of less than 0.5" cannot be explained by 
projection effects (Fig. 14 1. This implies that wider or larger flux 
contrast systems, if existing, cannot be individually disentan- 
gled from the pairs arising by chance alignment along the line 
of sight. 

Table |5] lists the 150 pairs separated by 5" or less and with 
/"bound ^ 0.99. The first column indicates the primary and sec- 
ondary IDs from Table [3] Columns. 2-3 and Cols. 4-5 list the 
H and magnitudes of the primary and secondary compo- 
nents, respectively. The separation of the pair is given in Col. 6. 
The two closest pairs detected have a separation of 0.24" and 
0.25" (only 600 AU at the Tr 14 distance) and AKs -0.6-0.7 
mag (Fig. 15 1. The closest probable companion to a massive star 
is a /Ts - 13.2 mag star at 0.4" from the Bl V star Trl4-19 
(pair ID #1530-1536 in Table |5]). Of the 31 stars brighter than 




Fig. 15. Close-up view in the Ks band image of the four closest 
pairs of Table 5] The displayed regions are 6"x6" and the sep- 
arations are al in the range of 0.24"-0.30". Their identifiers as 
given in Table 3] are overlaid on the images. 



Table 5. List of likely bound visual pairs (/"bound ^ 0.99). 







.prim 






di") 


656- 664 


17.490 


16.655 


18.313 


17.359 


0.24 


306- 299 


15.088 


14.798 


15.997 


15.380 


0.25 


1446-1451 


14.097 


13.323 


13.803 


13.560 


0.28 


902- 896 


17.697 


17.044 


18.709 


17.784 


0.29 


1024-1025 


14.842 


14.455 


15.941 


15.349 


0.29 


1091-1092 


16.165 


15.342 


15.979 


15.459 


0.29 


156- 157 


17.303 


16.715 


17.665 


17.084 


0.30 


1581-1579 


15.277 


14.653 


15.407 


15.030 


0.31 


1488-1491 


14.205 


13.941 


16.044 


15.530 


0.32 


1445-1453 


16.320 


16.014 


17.395 


16.981 


0.32 


1398-1389 


14.326 


13.931 


15.461 


14.942 


0.33 


1152-1149 


14.867 


14.498 


16.254 


15.229 


0.34 


1351-1343 


12.899 


12.542 


15.742 


14.876 


0.34 


828- 825 


17.984 


17.458 


19.244 


18.233 


0.34 


1780-1779 


16.527 


16.052 


16.829 


16.322 


0.34 


1320-1318 


15.386 


15.024 


16.267 


15.660 


0.34 


945- 955 


15.533 


15.122 


16.289 


15.724 


0.35 


1163-1155 


16.471 


15.241 


16.277 


15.674 


0.35 


1101-1104 


12.828 


12.320 


15.311 


14.589 


0.36 


1202-1200 


15.949 


15.449 


16.106 


15.506 


0.36 



note: The full version of the table is available in the electronic version 
of the paper or through the CDS; http://cds.u-strasbg.fr/ . 



/Ts = 1 1 mag, only six have high probability companions in the 
range of 0.4" to 2.5", i.e. 19%. Because of the observational 
biases discussed and of the filtering criteria applied, the statisti- 
cal distribution of the selected pairs is not representative of the 
underlying distribution and will not be discussed further 
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Fig. 16. Left: Average number of companions per star for a given 
maximum companion magnitude. Right: Cumulative distribu- 
tion functions (CDF) of the companion brightness. The central 
stars are taken in four ranges as indicated in the upper left-hand 
legend. Only pairs with separations in the range of 0.5"-2.5" 
have been considered. 



5.3. Companion frequency 

To search for variations in the companionship properties of dif- 
ferent sub-populations, we first computed the average number of 
companions per star, considering various magnitude ranges both 



for the central star and for the companions (Fig. 16 left panel). 
We applied the same colour and magnitude selections defined in 
Sect. |4.2| We further adopted an exclusion radius of 5" between 
each central star so that each companion is only assigned to one 
primary, preserving the independence of the different samples. 
In this particular case and to allow direct comparison between 
the different samples, we did not require that /Ts.prim < ^Ts.sec- 
For the central stars, we consider four ranges of magnitudes: 



n < 11, corresponding to M > 
11 < /*rs,prim < 13, correspond- 



- the massive stars : 6 < ^Ts.pr 
10 Mo MS stars); 

- the intermediate-mass stars : 
ing to 10 > M > 4 Mq stars; 

- the solar-mass PMS stars : 13 < /Ts.prim < 14, corresponding 
to 2.5 > M > 0.5 Mq PMS stars; 

- the low-mass PMS stars : 15 < /Ts.prim < 16, corresponding 
to 0.2 > M > 0.1 Mq PMS stars; 

where the isochrones of Fig. [7] have been adopted as guidelines 
for the mass conversion. The companion magnitude was chosen 
in the range of 9.5 < /Ts^sec < 18. Finally, we restrained our 
comparison to the 0.5"-2.5" separation regime. As we will show 
below (Sect. 6.1 1, the observational biases have a very limited 
impact on our results in that range. 

Under those assumptions, we found that massive MS stars 
have on average 3.8 + 0.5 companions, while solar-mass PMS 
stars have 3.0 ± 0.2 companions. The number of companions 
of intermediate-mass stars and of low-mass PMS stars are not 
significantly different from one another. Most of the difference 
is however found for /Ts.sec < 17: 3.5 ± 0.4 companions for MS 
stars against 2.5 ± 0.2 for lower mass stars. The difference is 



2.0 ■ 



1 .5 - 



1.0 ■ 



0.5 - 



0.0 [ 



6<Ks_prim< 1 1 
1 1 <Ks_prim<1 3 
1 3<Ks. 
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Fig. 17. Average number of companions per star as a function 
of the companion brightness. The central stars are taken in four 
ranges as indicated in the upper left-hand legend. A moving av- 
erage with a 2 mag bin has been used and the considered separa- 
tion range is 0.5"-2.5". The various envelopes give the Icr error 
bars. 



thus significant at the 2.5cr. This corresponds to a rejection of 
the null hypothesis that massive and lower-mass stars have the 
same number of companions with a significance level better than 
0.01. Because this is seen against the observational biases (see 
Sect. 16.2.2] ), this result is likely to be even more significant. 



5.4. Magnitude distribution 

We computed the distributions of the companion magnitudes 
for the various stellar populations considered above using a 2- 
mag wide moving average along the companion magnitude axis 
(Fig. 17 1. Most of the differences between the distribution func- 



tions of the massive stars and of the lower-mass stars result from 
the range /Ts.sec < 14. The more massive stars display thus about 
twice as many solar-mass companions as the lower mass stars. 
The situation is reverse for fainter, low-mass PMS companions, 
where the lower-mass stars tend to have more companions. The 
latter effect can however result from the difficulty to detect ex- 
tremely faint stars in the wings of the brightest stars. To allow for 
quantitative statistical testing, we also built the cumulative dis- 



16 



tribution functions (CDF) of the companion brightness (Fig 
right panel). Using a two-sided Kolmogorov-Smirnov (KS) test, 
one can reject at the 2cr level the null hypothesis that the high- 
mass and the solar-mass stars share the same companion CDF. 

5.5. Companion spatiai distribution 

To investigate the spatial distribution of the companions, we built 
the growth curves of the number of companions as a function 
of the separation. For the curves to be more robust against low 
number statistics, we concentrated on the total growth curve of 
a given stellar population rather than on the growth curve of in- 
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Fig. 18. Upper panel: growth curve for the companions of mas- 
sive stars. The plain line indicates the expected distribution for 
a uniform distribution of the companions in the field of view. 
Lower panel: same as upper panel for the companions of solar- 
mass stars. 




d_min (") 

Fig. 19. Cumulative distribution functions (CDFs) of the dis- 
tance from each artificial star to (i) the closest artificial star (plain 
curve), (ii) the closest source in field (dashed line), (iii) the clos- 
est bright star with /Ts < 1 1 (dashed-dotted line). 

estimate of the maximum contrast achieved as a function of the 
separation Sect. 5.1 The pairing properties are well described 
by chance alignment except for the closest pairs {d < 0.5 ") 



and for the pairs with similar magnitude components (Fig. 14 1, 
yet we identified 150 likely bound pairs (Sect. 5.2 1. Massive 
stars further tend to have more companions than lower-mass 
stars (Sect. 5.3 i. Those companions are brighter on average, thus 



more massive (Sect. 5.4 1. The spatial distribution of the compan- 



ions of massive stars is however not significantly different from 
those of PMS star companions (Sect. 5.5 1. The significance of 
those results is better than 2cr but no better than 3cr, and remains 
thus limited. The situation is however reminiscent of the case of 
NGC 661 1 where |Duche ne et aD(2001 ) found that massive stars 
are more likely to have bound companions compared to solar- 
mass stars. For Tr 14, the fov is definitely more crowded and 
probably more heavily contaminated, which seriously compli- 
cates both the companionship analysis and the interpretation of 
the results. 



dividual targets. We further limited the companions to masses 
above 0. 1 Mq, which roughly corresponds to a magnitude limit 
of - 16 mag and, as in the previous paragraph, we restrained 
our comparison to the 0.5-2.5" separation regime. 



Figure 18 compares the companion growth curves around 
high-mass and around solar-mass stars with the theoretical dis- 
tribution expected from random association with an underlying 
uniform distribution across the field. On the one hand, massive 
stars seem to have their companions statistically further away 
than expected from a uniform repartition. On the other hand, the 
growth curve of solar-mass PMS stars follows the expected trend 
from random association. Yet in both cases a KS test does not 
allow us to reject the null hypothesis that both realisations are 
compatible with the uniform distribution in the considered sepa- 
ration range. Similarly, a two-sided KS test does not allow us to 
claim that both growth curves are different from one another. 

5.6. Summary 

To summarize the results of this section, the closest pair of de- 
tected stars in our data is separated by 0.24" (~600 AU), in good 
agreement with the IQ of Sect.|2] Equation|6]gives an empirical 



6. Observational biases 

This section first describes an artificial star experiment designed 
to quantify the detection biases in the vicinity of the bright stars. 
It also presents two simple models that generalise the results of 
the artificial star experiment and provide an estimate of the im- 
pact of various observational biases on the results of this paper. 

6.1. Artificial star experiment 

Most of the difference in the companion properties of massive 
and lower-mass stars are found for companions in the range 
13 < ^Ts.sec < 14. In this section, we describe and analyse the 
results of an artificial star experiment that aims at better un- 
derstanding the limitation of our data in that range, justifying 
some of the choices made in the previous section. The PSF of 
the brightest sources indeed show strong and extended wings, 
which decrease the detection likelihood in the neighbourhood of 
a bright star As a by-product, the results of this experiment also 
provide an independent estimate of the photometric errors and 
of the completeness of our catalogue for the tested parameter 
range, but this is not our main purpose. 
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Fig. 20. Lower panel: Comparison between the input and output 
magnitudes for the artificial star experiment as a function of the 
separation the to closest bright star in the fov (/Ts < 11). The 
squares and diamonds respectively give Amag - /Ts^out - ^s,in 
and i/out - ^^in- Upper panel: same as lower panel for the colour 
of the artificial stars : A(H - K) ^ (H - Ks)om - (H - Ks)in. 
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Fig. 21. Artificial star recovery fraction as a function of the dis- 
tance d to a bright neighbour (Ks < 11). The plain line shows 
a linear interpolation between the detection and non-detection 
regime and is given by /i-ecov - {d- 0.35)/0.2, for d in the range 
0.35-0.55". 



The artificial star experiment follows a procedure similar to 
that presented in |Momany et al.| ( [2002|[2008| . In particular, stars 
with known H and magnitudes were simulated into the in- 
dividual H and images using the PSF of each image and 
taking into account the quadratic dependence of the PSF with 
the position in each frame. The entire reduction procedure was 
then repeated and the artificial stars were reduced as described 
in Sect. 121 

To the first order, stars brighter than /Ts = 1 1 in Tr 14 have 
masses of 10 Mq or more (see e.g. Fig. |7]i. As in Sect. |5] we 
adopted this limit for our massive star sample. We thus selected 
14 bright stars with 9 < A's < 11. Around each of them we 
simulated 50 companions with 13 < /Ts.in 14 and i/in - A's.in - 
0.35, spread in a 5" radius. The colour of the artificial stars were 
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Fig. 22. Detection probability as a function of the stellar magni- 
tude resulting from the effect of the crowding in the Tr 14 fov. 



chosen to reproduce the colour of typical PMS stars in Tr 14, 
which are the dominant type of sources in that magnitude range. 

Seven hundred artificial stars were thus simulated in the H 
and Ks images. To optimise the computation time of this com- 
plicated procedure, all the artificial stars were added simultane- 
ously. While this led to some heavier crowding than in the orig- 
inal field, this will be taken into account in the analysis. This 
further allows us to study the detection biases using artificially 
controlled pairs. 

6.1.1. Separation distribution 

Figure [19] shows the cumulative distribution of the separations 
of the artificial stars with respect to one another, with respect 
to the closest source in the field and with respect to the closest 
bright source in the field. It shows that we are able to investigate 
various ranges of separation, from the crowding of the artificial 
stars among themselves to the effect of field density. We first 
focus on the closest detections. The main results of the artificial 
star experiment in this respect are : 

- The closest recovered pair of artificial stars is separated by 
t^sep = 0.2", in good agreement with the IQ derived earlier, 

- The recovery fraction of artificial pairs for which both com- 
ponents are further away than 0.6" from any source in the 
image is better than 0.99 for <isep > 0.24", 

- Excluding all pairs of artificial stars with t/sep < 0.3", the 
closest separation between a recovered artificial stars and a 
star at least as bright is 0.23", 

- Similarly, the closest separation between a recovered artifi- 
cial star and a massive star (K^ < 11) is 0.38", in perfect 
agreement with what we found in our data (Sect. 5.2 1, 

- All in all, there is a good agreement between the results 
of the artificial star experiment and the maximum reachable 
contrast as a function of the separation described empirically 
by Eq.|6] 



As a first result, the artificial star experiment allows us to vali- 
date up to AKs - 5 mag at least the empirical contrast vs. sepa- 
ration function introduced earlier This result will be used below 
to build a first-order analytical model of the observational biases 
affecting our data. 
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6.1 .2. Recovery fraction 

To estimate the recovery fraction of solar-mass PMS stars in the 
wings of the brighter massive stars, we first excluded all the close 
artificial pairs from our analysis. We also exclude all the artificial 
stars that fall closer than 0.3" from any source fainter than Ks - 
11 in our data. With this we eliminate the uncertainties due to (i) 
crowding among the artificial stars and (ii) confusion between 
the artificial stars and the numerous field sources in our data. 
Almost 500 artificial stars are left, providing a decent coverage 
of the parameter space. 



Figure 21 shows the recovery fraction as a function of the 
separation to the closest bright star. As mentioned earlier, the 
first detection occurs for a separation of 0.38". For t/sep > 0.53", 
all the artificial stars are recovered. Between 0.38" and 0.53", 
the recovery fraction is approximately 0.4 and remains constant 
over the interval. 



6.1.3. Photometric uncertainties 

The comparison between the artificial star input (/Tsjn) and re- 
covered (A's,oiit) magnitudes provides us with a more realistic es- 
timate of the photometric errors (Fig. 20 1. It reveals that the ac- 
curacy of both the retrieved magnitudes and of the colour term is 
significantly affected within x I" from a bright star, the fainter 
companion being up to 0. 1 mag too red. It also reveals a slight 
systematic shift in the retrieved colour and magnitudes, even at 
larger distance. Rejecting the few significantly deviant points re- 
sulting from the crowding, we obtained for the artificial stars 
more distant than 1.5" from the closest bright star 



(H-Ks)o 



Horn - Hm = 0.008 + 0.020 
^^s.out - ^s,in = 0.008 + 0.026 
-(H-Ks)in = -0.021+0.019, 



(10) 
(11) 
(12) 



Those are slightly larger deviations that the formal errors of the 
PSF photometry (Fig. |6]l. The systematic increase of the H- 
colour when getting closer to the bright companion could result 
from a slightly better IQ obtained in the Ks band compared to 
the H band, and thus a better subtraction of the bright star wings 
for a given separation. 

6.2. Analytical models 
6.2.1 . Impact of crowding 

Given the very steep transition in the detection probability once 
the separation to a brighter source increases, and because the 
artificial star experiment generally agrees with the empirical de- 
tection limit given by Eq.|9] one can develop a very simple model 
to estimate the impact of the crowding in the field. It relies on the 
following hypotheses : 

- A brighter star is always detected if falling on top of or very 
close to a fainter one, 

- Non-detection is the consequence of shadowing by brighter 
stars in the fov, 

- The spatial distribution of the stars in the field is random. 

The detection probability of a source of a given magnitude 
can then be written as 



^detect, (^s) ~ 1 



(13) 
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Fig. 23. Detection probability model for pairs with separation in 
the range of 0.5"-2.5" and for various brightness ranges of the 
central star. 



where d(AKs) is given by Eq.[9]and where A^v is the area of the 
considered fov. In Eq. [T3] th e sum is performed over all sources 
brighter than /T". Figure : 



22 



, shows the resulting detection proba- 
bility. According to this model, completeness level of 0.99, 0.95 
and 0.90 are reached for stars brighter than Ks = 14, 16 and 
17 respectively. While our model is certainly too crude to ade- 
quately describe the faint end, it still shows that the crowding 
has a very limited impact in Tr 14. In particular, this indicates 
that the star counts used to fit the Tr 14 profile and to estimate 
the cluster core mass are unlikely to be significantly biased. 

6.2.2. Companion detection tliresliold 



In Sect. |5.3| to 5.5 we focused our analysis to the 0.5-2.5" sepa- 
rations. In this section, we develop again a simple model to better 
quantify the impact of the observational biases in that range. As 
above, we use Eq. |9]to define the region where a star outshines 
close fainter neighbours given the magnitude contrast of the pair. 
In particular, the detection probability of a neighbour of magni- 
tude Ks^sec in the vicinity of a /rs,p,-im central star can be mod- 
elled as the ratio between the area where one of the components 
does not outshine the other and the total area considered. For an 
annulus region /?min-^max around the central star, the detection 
probability can thus be written as 



r 1.0, 

1.0- 
0.0, 



ifd< R" 



Rm 



min ^ d < ^max' 



ifR 



(14) 



Figure 23 compares the detection probability in the 0.5-2.5" sep- 
aration range for the various central star-brightness ranges used 
earlier For the bright primary interval, the curve displayed re- 
sults from the average of Pdetect. computed individually for all 
the bright primaries considered to build Figs. 16 and 17 The 



ifov 



central magnitude of the primary intervals has been used for the 
other three categories. 

Our results show that the detection of the companions to 
intermediate-mass stars (11 < Ks < 13) and to solar-mass PMS 
stars (13 < Ks < 14) stars is mostly unaffected. As expected, the 
lower mass PMS (15 < ^Ts.prim < 16) are less likely to be found 
close to a bright stars. However, taking into account the number 
of bright stars in the field, one can prove this to be completely 
negligible. The largest biases are affecting the companions of 
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Fig. 24. Evolution of Amsr with the adopted magnitude limit 
^s.min for the massive star sample. The circles and triangles in- 
dicate the results obtained with and without the colour and mag- 
nitude selection criteria of Fig. |8] The upper and lower values 
give the number of stars considered in each step and for both 
samples. The dotted line indicates Amsr - 1> no mass seg- 
regation. 



bright stars. Pdetect. is passing below 0.9 at Ks - 14 and below 
0.5 at Ks ~ 17. One can further show that up to one companion 
per star is likely lost at /Ts.sec < 16 in Fig. 16 



7. Mass segregation 

As introduced in Sect. |4.3| the more massive MS stars seem more 
concentrated towards the cluster centre than the lower mass PMS 
stars. The best-fit EFF87 profiles (Table 4]i confirm that this dif- 
ference is indeed significant at the 4cr evel. This could be in- 



terpreted as a hint for mass segregation, although Ascenso et al. 
( |200 9 ) warned against hasty conclusions because numerous ob- 
servational biases are actually favouring the detection of mass 
segregation, even in non-segregated clusters. 

[Allison et al.| ( |2009| l recently introduced an alternative 
method to investigate mass segregation, which is insensitive to 
biases like the exact location of the cluster centre, and less sen- 
sitive (although quantification is still lacking) to the incomplete- 
ness effects. Their method compared the minimum spanning tree 
(MST), the shortest open path connecting all points of a sample, 
of the massive stars to the equivalent path of low mass stars (see 
Allison et al.| for a full description of the algorithm). The mass 
segregation ratio (Amsr), i-C- the ratio between the average ran- 
dom path length and that of the massive stars, allows them to 
quantify the deviation between the massive star sample and the 
reference sample. Following their approach, Fig.|24]displays the 
evolution of Amsr with the adopted magnitude limit (/Ts^min) for 
the massive star sample. The reference distribution consists of 
stars with 14 < Ks < 16 and was drawn 500 times from our cat- 
alogue, with each sample containing the same number of stars as 
found in the bright sample. The dispersion obtained gives us the 
error bars on Amsr as displayed in Fig. 24 In the above proce- 
dure, we deliberately remained far from the limiting magnitude 



of our catalogue to minimize the completion biases. The method 
is in principle still aff'ected by crowding and by the shadowing 
in the vicinity of bright stars. We showed in Sect. 6.1.2 however 



that the former effect had a very limited impact down to A's < 16 
at least. The effect of the shadowing in the vicinity of bright stars 
is more difficult to estimate, although one can expect that the ab- 
solute number of 14 < Ks < 16 stars lost is proportionally very 
small compared to the number of stars in that interval. This re- 
sults from the low number of bright stars and from the limited 
radius at which they can outshine a fainter neighbour. 

We computed the MST first with and without the colour se- 
lection defined in Sect. 4.2 (Fig. 24 1. Focussing on the most 
probable members (thus applying the colour selection), we 
found some indication of mass segregation down to ATs^min ~ 
10.5 mag at the 1.5cr level. As expected, the degree of mass seg- 
regation is increasing with the average brightness of a sample, 
thus with the mean stellar mass. The largest mass segregation 
ratio is obtained for the few brightest, most massive stars after 
applying the membership selection. Yet the 1.5cr confidence of 
this result remains at the limit of the detection. 

Because the sensitivity of the MST to the completeness of a 
sample is not fully understood (Allison et al. 2009), we cannot 
draw firm conclusions. We note however that two independent 
methods, profile fitting and MST analysis, both point towards 
mass segregation, as can be expected for the most massive stars 
in such a cluster. 

Obviously our observations only allow us to investigate the 
current mass segregation status of the cluster and we cannot dis- 
tinguish whether this segregation, if confirmed, is primordial or 
is the product of early dynamical evolution. Given the expected 
cluster mass and size and its stellar contents as obtained in 



Sect. 4.3 after colour selection, we estimated the typical dynami- 
cal friction time-scale fdf of 10 Mq and 20 Mq stars (correspond- 
ing to resp. /Ts ~ 1 1 mag and 9.5 mag in our data). Following 
Spitzer & H'Etl ( [T971| and Porte gies Zwart & McMm an (2002|l 
, we obtained fdf ~ 7.2 x 10^ yr and 3.6 X 10^ yr respectively. 
Considering the estimated age of the cluster, 3-5x10^ yr, the dy 
namical friction time-scales agree with the results of Fig. 24 



where mass segregation begins to appear somewhere between 
10 Mgand 20 Mq. As a consequence, if mass segregation is con- 
firmed, it does not need to be primordial but can probably be 
explained by dynamical evolution. 



8. Summary and conclusions 

Using the ESO MCAO demonstrator MAD, we have acquired 
deep H and Ks photometry of a 2' region around the central part 
of Tr 14. The average IQ of our campaign is about 0.2" and the 
dynamic range is about 10 mag. The image presented in Fig. [T] 
is by far the largest AO-corrected mosaic ever acquired. 

Using PSF photometry, we investigated the sensitivity of 
faint companions detected in the vicinity of bright sources. We 
derived several empirical relations that can be used as input for 
instrumental simulations, to estimate the performance of AO 
techniques versus seeing-limited techniques or, as done later in 
this paper, to build first-order analytical models of the impact of 
some observational biases. In particular, the contrast vs. sepa- 
ration limit has been validated over a 5 magnitude range by an 
artificial star experiment. 

Despite a probably significant contamination by field stars, 
the Tr 14 CMD shows a very clear PMS population. Its location 
in the CMD can be reproduced by PMS isochrones with con- 
traction ages of 3 to 5 X 10^ yr. Interestingly, Tr 14 cannot be 
significantly further away than the distance obtained by |Carraro| 
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|et al.| ( |2004| l i.e., 2.5 kpc, as this would result in an even ear- 
lier contraction age. We derive the surface density profile of the 
cluster core and of different subpopulations. For stars brighter 
than /Ts = 18 mag, the surface density profiles are well re- 
produced by EFF87 profiles over our full fov, and we provide 
quantitative constraints on the spatial extent of the cluster and 
on its stellar contents. Adopting the core-halo description sug- 
gested by Ascenso et al. (2007 ), we report that the transition be- 
tween the core and the halo is not covered by our data, implying 
that the core is strictly dominating the density profile in a radius 
of 0.9 pc at least. Using colour criteria to select the most likely 
cluster members, the density profiles of the more massive MS 
stars are best described by a power-law (or, equivalently, by an 
EFF87 profile with a very small core radius). 

We also investigated the companionship properties in Tr 14. 
We showed that the number of companions and the pair asso- 
ciation process is on average well reproduced by chance align- 
ment from a uniform population randomly distributed across the 
field. Only stars with a brightness ratio close to unity or with a 
separation of less less than 0.5" cannot be explained by spuri- 
ous alignment and are thus true binary candidates. This does not 
imply that large light-ratio and/or wider pairs do not exist, but 
rather that they cannot be individually disentangled with statisti- 
cal arguments. Still, 19% of our massive star sample have a high 
probability physical companion. 

Focusing on the 0.5"-2.5" separation range, where the obser- 
vational biases are unable to invalidate our results, we compared 
the companion distributions of massive stars with those of lower 
mass stars. In Tr 14, the high-mass stars (M > 10 Mq) tend to 
have more solar-mass companions than lower-mass comparison 
samples. Those companions are brighter on average, thus more 
massive. Finally, no difference could be found in the spatial dis- 
tribution of the companions of low and high-mass stars. 

Lastly, we employed the MST technique of [Alhson et al. 
( |200 9 ) to investigate possible mass segregation in Tr 14. Again 
we found marginally significant results (at the 1.5cr level), sug- 
gesting some degree of mass segregation for the more massive 
stars of the cluster (M > 10 Mg) . Although the sensitivity of the 
method to incompleteness is still not fully quantified, we note 
that early dynamical evolution can reproduce the observed hints 
of mass segregation in Tr 14, despite the cluster's young age. 
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